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ABSTRACT 


This report examines the problem of calculating solar panel power output 
at launch and during a space mission. It also examines the major sources of 
uncertainty and error in predicting the post launch electrical performance of 
the panel. A general discussion of error analysis is given. Examples of 
uncertainty calculations are included, A general method of calculating the 
effect on the panel of various degrading environments is presented, with 
references supplied for specific methods. A technique for sizing a solar panel 
for a required mission power profile is developed. 
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SECTION I 


INTRODUCTION 

This memorandum presents calculational procedures for predicting the 
power output of a solar panel throughout its life. It shows how to determine 
the power output of a freshly manufactured solar panel, of a panel immediately 
after launch, and of a panel after exposure to the rather poorly known space 
environments, 'which may include radiation, ultraviolet exposure, temperature 
excursions and changes in incident solar intensity. Uncertainties in making 
the predictions are detailed and a method of calculating the overall uncertainty 
in panel performance is presented. The prediction uncertainty is the main 
subject treated in this paper. The problem is separated into two broad cate- 
gories: (1) prediction of panel performance immediately after launch and 
(2) prediction thereafter. 

Post launch environments affecting the panel output which are amenable 
to treatment by the techniques in this paper include solar intensity, panel 
temperature, electron and proton irradiation arising either from the Van Allen 
belts or solar flares, ultraviolet degradation of the cover slide optical system, 
micrometeoroid irradiation, and the effect of thermal cycling. 

A summary of the factors treated in this paper which contribute to the 
uncertainty in predicting array space performance include: 

1) Errors in cell measurement which contribute to cell matching loss 

2) Effects of cover slide filtering 

3) Manufacturing losses due to soldering, wiring, etc. 

4) Calibration of the balloon flight standard solar cells 

5) Prelaunch Table Mountain (Earth sunlight) measurements of space 
solar arrays 

6) Uncertainties in the JPE computer program for predicting panel 
performance in space 
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7) The following damaging environments to which the array will be 

exposed in space. (The effect of these environments is a function 

of the spacecraft mission parameters and the solar array design): 

a) Ultraviolet radiation 

b) Micrometeoroid impact 

c) Van Allen belt radiation 

d) Solar flare protons 

e) Thermal cycling 

All the uncertainties due to prelaunch measurements and calibration 
are combined to calculate an uncertainty in the predicted panel output just after 
launch. The methods of predicting solar panel output as panel temperature 
and incident solar intensity change used in the JPL computer program M132000 
are described and criticized. The errors in using the methods in the program 
are analyzed, tabulated and combined. An example of an error calculation for a 
panel like the Mariner '71 Mars orbiter solar panel is presented. Recommend- 
ations are given which are designed to decrease the panel assembly losses and 
reduce the panel prediction uncertainties. 

Using the JPL Ml 32 000 computer program to predict the panel output, an 
uncertainty in the prediction of I , V , and P may be computed as 

SC O C H13-X 

follows; 



'dv \ 2 / . av \ 2 / . av \ 2 


oc 
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dV V \ Z /dl I \ 2 

oc oc \ , / sc sc \ 

V V ) + \1 1 ) 

oc mp / \ sc mp / 


The terms appearing in the above equations are defined below: 

I = cell current during parametric data acquisition at (H, T ), 
mA 

I = panel current, mA 

P 

H = solar intensity at which solar cell output is predicted, 
mW / cm^ 

dH = solar intensity uncertainty including simulator intensity 
uncertainty when the parametric solar cell data used in 
Ml 32 000 was acquired and the uncertainty of the solar 
intensity incident on the spacecraft. 

n = temperature coefficient of short circuit current at the 

nearest tabulated (H, T) point in M132000 to the desired 
intensity and panel temperature, mA/°C. 

a - same as a^ t but at the desired intensity and panel 
temperature 

dT^_ = cell temperature uncertainty during parametric data 
acquisition °C 

dT_^ = panel temperature uncertainty, °C 

n = number of cells in parallel on the panel 

/dl\ 

lyj = uncertainty in measuring solar panel current at Table 

™ Mountain 

Sj = statistical 95 percent confidence limit of the parametric 
current data, mA 

= current prediction error of M132000, mA 


JPL Technical Memorandum 33-673 



V ocp = panel open circuit voltage, mV 


V QC ~ single cell open circuit voltage, mV 


m = number of cells in series on the panel 


- uncertainty in measuring the solar panel voltage at 
TM Table Mountain 


S = statistical 95 percent confidence limit of the parametric 
v 

voltage data, mV 


e = voltage prediction error of Ml 32 000, mV 


^m. = P ane ^ maximum power, mf 


V mp - voltage at maximum power, mV 

I = current at maximum power, mA. 
mp 

In the example calculations the uncertainties in predicting I gc , V qc , and 

P at 100 mW/cm 2 and 0°C were calculated for a Mariner '71 panel. At 
max 

the 95 percent confidence level they are 2. 49 percent for I gc » 2. 26 percent 

for V and 3. 82 percent in 

oc max 
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SECTION II 


PANEL ASSEMBLY POWER LOSSES 

Solar cells for a particular project are purchased to a specification which 
imposes some restriction on the power output of each cell. The restriction is 
usually applied to solar cell current at a constant voltage because this is a 
relatively easy measurement to make on the assembly line and can be a fairly 
good approximation to the cell maximum power. A distribution function for 
the cell power can be obtained from a proper sampling of the total population 
of solar cells manufactured for that project. Means and standard deviations 
can be estimated. Such a distribution for Lots IB through 5 of the Mariner 9 
(interchangeably referred to in this paper as Mariner 1 7 1 ) solar cells is shown 
in Figure 1. These cells were 1 to 3 ohm- cm and measured 2 x 2 x 0. 046 cm. 

The distribution shown is for 1000 cells sampled at the manufacturer’s plant 
and is plotted for cell current at a voltage of 0.485 volts. 

At this point several things are done to the solar cells in the process of 
assembling them into a panel which change the power distribution function. Data 
regarding the changes to the distribution function as the cells are processed is 
difficult to find. The cells are sorted by the manufacturer into output current 
groupings of 2 mA. A sampling of the output lot is taken at the manufacturer's 
facility for the purpose of establishing the overall current distribution, averages, 
standard deviations, and to determine if the cells will pass their required environ- 
mental tests. If a lot of cells is satisfactory, the assembly of submodules is 
begun. A submodule consists of a number, p, of cells (4 or 5 on Mariner ! 71) 
which are connected in parallel. The cells comprising a submodule are selected 
so that the submodule current output (at reference voltage) will be p times the 
average output current. If each submodule has nearly the same current output 
near operating voltage, no one single submodule will act as a load when the sub- 
modules are later connected into a series string. Coverglasses are now applied. 

In order to insure that the submodules will have matching current outputs even 
after the application of the coverglasses, the current output of each sub- 
module is remeasured. The practice at this laboratory has been to measure 
the submodules in tungsten light after glassing. The submodules are then 
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connected by a serie s /parallel arrangement into sections and the sections 
connected together to match the input impedance requirements of the spacecraft 
power conditioning system. 

It is clear that several or all of the above steps change the current output 
of the solar cells. The maximum power output of a typical solar panel built 
for JPL missions is known to be some 6 percent lower than that computed by 
multiplying the number of cells on the panel by the average power output of a 
single cell. Some data on the losses due to glassing will be presented in the 
following paragraphs, but no data could be found on the effect of soldering inter- 
connects onto the cells, the effect of the cell loading due to the mismatch in 
submodule current output or in the errors made by matching submodules after 
glassing on the basis of performance in tungsten light. 

A large portion of the manufacturing loss occurs when the coverglasses 

are attached. Cell vendors were cooperative in furnishing data regarding 

these losses. One vendor furnished complete I-V curves for 128 solar cells 

(N/P, 7-14 ohm-cm, 2 x 2 x 0. 025 cm (10 mils thick)) both before and after 

installation of 0. 15 mm (6 mil) microsheet coverglasses; the other vendor 

furnished data on 100 solar cells (N/P, 1-3 ohm-cm, 2 x 4 x 0. 03 cm (12 mils 

thick)) before and after glassing with 0. 15 mm (6 mil) 7940 fused silica cover- 

glasses. All cells were coated with a SiO coating. The short circuit current 

distribution for both types of cells both before and after glassing is shown in 

Figures 2 and 3. Superimposed on the histograms of current distribution are 

Gaussian curves having the same means and standard deviations. Chi-square 

tests at the 95 percent confidence level of the 2x2 cm cell data showed that 

the current distributions of the bare cells fit a Gaussian distribution, but did not 

after glassing. The 2x4 cm cells fit Gaussians both before and after glassing. 

The average I gc loss for the 2 x 2 cells is 3. 658 mA or 2. 66 percent when 

glassed with 6 mil microsheet. The average I loss for the 2x4 cells is 

s c 

5. 577 mA or 1. 97 percent. The percentage loss for cells covered with 7940 
fused silica is less, as is to be expected, because of the lesser light transmission 
of the microsheet. Figure 4 is a scatter plot of the 2 x 2 cm cell short circuit 
current loss as a function of initial short circuit current. The losses have a 
wide scatter and are unpredictable. The standard deviation in cell loss is 
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1.28 mA. Figure 5 is a similar plot for the 2x4 cm cells. Figure 5 is 
plotted with both scales twice that of the scales used in Figure 4 for the 
2x2 cells. Here the scatter is less. The standard deviation of the cell loss 
is 1.43 mA. Both vendors offered the opinion that the scatter in the loss data 
is primarily due to variations in the SiO coating which ailect the absorption of 
light in the coating as well as its index of refraction. There are too many uncon 
trolled variables in these data to draw conclusions as to whether cells glassed 
with 7940 have less current loss scatter than cells glassed with microsheet. 

The data on the 2x4 cm cells was obtained under unusually intense scrutiny of 
the people applying the co ver glass e s , while the same may not have been true 
for the 2x2 cm cells. The cells are of different resistivities, different sizes, 
may have had different cutoff frequencies for the ultraviolet filters, different 
coverglass adhesives, or at least different batches of the same adhesive, etc. 

It is clear however that there is a large scatter in the current loss data and 
cells which were matched on the basis of preglassing current output for the 
purpose of assembly into modules will have a significant mismatch after 
glassing. This loss and the resulting mismatch will account for at least 
1.97 percent of the manufacturing losses when 6 mil coverglasses are used 
and more when thicker coverslides are used. It is possible that this trans- 
mission loss and mismatch thereafter may account for the entire 6 percent 
loss seen in recent Mariner panels when 20 mil coverslides were used. 

It is recommended that the cells be glassed prior to measuring the output 
current for matching purposes. One of the vendors also suggested that the 
use of a TiO x anti reflection coating would cut down the size and scatter of the 
current loss since this material has an index of refraction which more nearly 
matches Si than does SiO. It is therefore further recommended that an investi- 
gation be initiated for the purpose of determining why there is so much scatter 
in the current loss due to glassing, and to try to find ways, including the use 
of TiO x anti reflection coatings, to reduce the scatter. 
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SECTION in 


GENERAL ERROR DISCUSSION 

A large portion of the remainder of this paper will involve discussions of 
measurement errors and their propagation. This section will discuss philosophy 
and methodology to be used in that regard. 

An excellent discussion of errors and uncertainties is given by Kline and 
McClintock in Reference 1. A brief summary of a few of their points follows. 

Uncertainty is defined to be the possible value the error might have. 

Note the distinction between error and uncertainty. The difference between the 
true value in an observed variable and the observed value is the error in that 
observation. The uncertainty, however, is what an observer thinks the error 
might be. 

An error frequency distribution function can be constructed to describe 
the errors if a large number of measurements of a variable can be made. Such 
error frequency functions are often non- Gaussian, but they almost always show 
that small errors are more likely than big ones, positive and negative errors 
are about equally probable, and there is no finite upper limit to the possible 
size of an error . In spite of this, design documents often specify maximum 
errors or uncertainties. Such numbers should be regarded as highly suspect. 

It is seldom possible to measure the uncertainty distributions of all the com- 
ponents contributing to the uncertainty of a variable, and another method of 
describing the uncertainty is required. One method of quoting the best estimate 
of the value of a variable and the associated measurement error is to give the 
mean m (arithmetic mean of observed values) and an uncertainty interval w 
based on specified odds, b. Thus a variable might be quoted as 

m ± w (b to 1 ) 

which means that the observer is willing to wager odds of b to 1 that the true 
value lies between m - w and m + w. The larger the w, the longer the odds one 
can safely wager. Expressed another way, this interval of width 2w centered 
on m has a probability of b/(b 4- 1) of containing the true value of the variable. 
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Let an observation R be the function of n independent variables 


R - R(v^, v 2 , . , . v n ) (1) 

If the variables v. are normally distributed, with uncertainties w. associated 
1 1 

with each v., and all w. have the same odds, then the uncertainty interval w^ 

1 i R 

for the result having the same odds is given by 



This equation gives results that are very nearly correct even when used with 
uncertainty distribution functions that are grossly non-normal (e. g. , 15 per- 
cent error for a traingular distribution). This is quite a remarkable result in 
that it holds for distributions which have uncertainty distributions with finite 
limits such as the triangular distribution. In view of the fact that uncertainty 
intervals for the variables are not usually known to better than 50 percent, the 
use of equation (2) is not unreasonable. 

To make use of the uncertainty interval based on odds as discussed by 
Kline and McClintock, we make use of the statistical confidence limit. Suppose 
we have a sample population distributed normally, whose true mean p and true 
variance c r are somehow known. If we draw a sample of size n at random 
from this population and calculate a mean x from this sample, we can calculate 
a confidence interval which has a probability (1 -a) of containing the true 
sample mean. This interval is commonly called a 100(1 -cp) percent confidence 
interval. The bounds of the interval are given by 

s±z f (t) < 3) 

In gambling language, the odds are b to 1 that this interval contains the true 
mean p. b and n are related by b = (1 - a) /a. The a used here is a common 
notation symbol in statistics texts to describe confidence limits and is not to be 
confused^with the a commonly used by photovoltaic s engineers in describing 
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the temperature coefficient of solar cell short circuit current. 

value in the standard normal distribution such that the probability of a random 

deviation greater than ± Z a jz a * i* e * 

CO 

f(Z) dZ = 

2 



where 


f(Z) = 


\/2tt 



and Z 


x -K 
cr 


If we are using 95 percent confidence limits, oe - 0. 05 and Z = ^ 96, 
and our interval is X ± 1. 96c r/\/iu This interval has odds of 19 to 1 for con- 
taining the real population mean p. 

Suppose our samples are drawn from a population just like the one in the 
example above except that the variance is not known . In this event, which is 
the usual situation encountered in the real world, we typically estimate the 
variance with a value s found from a random sample of size n as follows; 


2 

s 


1 

n-1 



i - 1 


or equivalently 



(4) 


It is most important to note that the factor of 1. 96 used in the previous example 
will no longer insure 95 percent confidence limits. A bigger factor is needed 
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2 • 2 

because s is only an approximation for cr . In this case the 100(1 -o') percent 
confidence limits are 


X ± 



( 5 ) 


Here, ^ctj2, n-1 integral of the Student t distribution (from t to co) for 

n -1 degrees of freedom. This distribution is tabulated in most statistics 
books. For n - co, t q ^5 “ coinciding with our previous result for 

known variance, and suggesting that if we could draw an infinite sample size we 
^ould compute an accurate s. If we drew a sample of size 5j however, we would 
have t ^ = 2. 776 as the multiplier of s/\/n7 

One often hears the terms lo~, 2cr, or 3a* levels or limits used in a rather 
cavalier fashion in conjunction with engineering data. The implication usually 
is taken that formula (3) can be used to set confidence limits. This is rarely 
the case however, and equation (5) must be applied. Further discussion of 
these points along with examples of their application can be found in 
References (2) and (3). 
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SECTION IV 


CELL MEASUREMENT UNCERTAINTY AT THE CELL 
MANUFACTURER'S FACILITY 

We consider that the solar cell manufacturer has set up his simulator to 
measure the current output of a batch of solar cells from his assembly line. 

He is following the requirements of the specifications such as those detailed in 
Reference 4. This is done by using a standard solar cell called a secondary 
working standard which is placed in a temperature controlled holding fixture 
with integral spring loaded contacts. The standard cell is brought to the 
required temperature of 28°C and the solar simulator intensity adjusted until 
the short circuit current output of the standard cell attains its calibration value 
for one solar constant. But, due to uncertainties in the calibration and measure- 
ment process, he has probably not actually achieved the desired simulator 
intensity, H^, or the desired temperature, T^. Instead the simulator is set 
at an intensity of H = H q + dH and the standard cell is set at a temperature of 

T , = T + dT. If we knew absolutely what H, and T, were, we could correct 
I o 1 1 

the simulator intensity by changing the short circuit current of the standard 
cell an amount 


dl = 


31 

3H 


dH + 


31 

3T 


dT 


or 


dl dH , dT 

T = h~ + o T" 


( 6 ) 


if dH and dT are small. Equation (6) assumes that I = const, x H, so that 

(l/I) (3I/3H) = 1 /H, a is the temperature coefficient of short circuit current 

at intensity H and temperature T . (Uncertainties in the value of a contrib- 
o^o o 

ute negligibly to dl/I and are neglected. ) An alternative interpretation to equa- 
tion (6) is that the uncertainty in intensity may be written as 


jtj. j T a dT 
dH _ dl o 

H I " I 
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with dl representing the possible error the standard cell current may have 
because of uncertainties in the chain of measurements involved in its 
calibration. 

Uncertainties in the current calibration begin to arise with the calibration 
of the balloon flight standard cells (BFS) which are flown by JPL on a balloon 
at altitudes of 24. 38 m (80, 000 ft) or higher. These cells are loaded with integral 
precision 1 ohm resistors (0. 1 percent) and the voltage drop telemetered to a 
receiving station. One resistor is used on a 1 x 2 cm cell and 2 resistors in 
parallel are used on the 2x2 cm cells. The cells are recovered after the 
balloon flight and used in the laboratory as a transfer standard. It does not 
matter what the accepted value of the solar constant at air mass zero (AMO) 
is, nor does it matter what the actual value of the load resistor is as long as 
it is stable, has a low temperature coefficient, and loads the cell near short 
circuit condition. If the BFS cell is put into a solar simulator beam whose 
spectrum is a reasonable approximation to that of the sun, and the .simulator 
adjusted to give the same voltage drop across the BFS load resistor as it had 
at 24. 38 (80, 000 ft. ) with the cell at the same temperature, we have. used the 
BFS to transfer the AMO intensity into the laboratory. In performing this 
transfer, the sources of error are errors in the voltage read-out and tempera- 
ture control system during the ballon flight, and voltmeter and cell temperature 
measurement errors at time of use in the laboratory. 

The solar cell vendor does not use a BFS cell to set the intensity of his 
simulator, however. The JPL specification of reference 4 allows the vendor 
to use a cell called a secondary working standard (SWS) to set the simulator 
intensity. This SWS cell is calibrated against a BFS cell, but in doing so, a 
certain loss of accuracy occurs. A freshly calibrated BFS cell is considered 
to be accurate to ± 1/2 percent (95 percent confidence limit), but degrades as 
it is used extensively in the laboratory to an estimated ± 1 percent. The SWS 
calibration is taken to be ± 2 percent (95 percent confidence limit). The 
temperatures of the BFS and SWS cells are commonly held to ± 1°C which can 
introduce approximately 0. 1 percent uncertainty in the intensity when a 2 ohm 
cm standard cell is used. The BFS and SWS cells are read with a digital 
voltmeter assumed to be accurate to ± 0. 1 mV out of approximately 6 0 mV. 

An additional uncertainty of 0. 2 mV is assigned due to difficulties in reading 
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the cell current because the simulator intensity has short term fluctuations 
caused by the dancing of the arc in the xenon arc lamp. We therefore assign 
an uncertainty of 0, 5 percent to the voltmeter readings. Using a SWS cell to 
set up the intensity of the simulator, and if the SWS cell reads exactly the 
calibrated value, a total uncertainty in the simulator intensity may be calcu- 
lated using equation (2) as follows: 


dH 

H 



t (0. 001) + (0.005) 


n 


1/2 


X 100 = 2. 064% 


( 8 ) 


In addition to the BFS and SWS calibration uncertainties, the solar cell 
manufacturer is allowed to set his simulator intensity to within ± 2 percent of 
AMO (Reference 4).‘ It is not known how often the manufacturer recalibrated 
the light source intensity during a day's series of measurements, but the 
manufacturer's specification for the X-25 solar simulator quotes a ± 1 percent 
intensity variation to be possible over a time period of 6 0 minutes. Assuming 
these additional uncertainties to be 95 percent confidence limits, we can 
compute via equation (2) the total uncertainty in the simulator intensity to be 


dH = 
H 


[c°. 


02064/ + (0. 02) + (0. 01) 


2 ] 


1/2 


X 100 = 3. 04% 


(9) 


The solar cell vendor now proceeds to read the current output of his batch 
of solar cells. The uncertainty in the cell current output is given by equation (6), 
with the additional uncertainties of the current measuring device factored in. 

The current is read by using a digital voltmeter to read the voltage drop across 
a precision resistor, R^. The current uncertainty introduced by the voltmeter 
and resistor is 


dl 


m 


m 


,,, dR 
dv c 


V 


R 


( 10 ) 


where dl^ represents the uncertainty inherent in the current meter alone. We 
assume the manufacturer uses a 0. 1 percent resistor and a digital voltmeter 
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with a calibration uncertainty of ± 1 digit in the least significant place (1 part 
out of 1200 when measuring currents near 120 mA at the 485 mV point). Again 
the xenon arc light will flicker and give an uncertainty of approximately 
2 digits in the least significant place. The total uncertainty of the current 
reading device is computed by the scheme of equation (2) to be 



m 


( 0 . 001 ) + 



0. 27% 


(ID 


Using equations (2), (6) and (11), the uncertainty in the solar cell current at 
V = 485 mV is computed as follows: 


dl 

I 


G 


(0. 0304) + 


( o.llx l V 

\ 120 ) 


-, 1/2 


+ (0. 0027)' 


X 100 = 3. 05% 


( 12 ) 


assuming that the manufacturer is able to maintain the cell temperatures 
within 1°C. The nature of this uncertainty is to define a 95 percent confidence 
band on the solar cell I-V curve by adding ± Al to each point on the curve such 
that AI = 0. 0305 x I . Such a band is shown on a typical Mariner ! 71 cell in 
Figure 6. 

The spectral content of the test light source can also affect the output of 
a solar cell. However, according to tests reported by R. Mueller {Ref. 5), 
quite wide variations in simulator spectral content did not induce directly 
dependent changes in solar cell output when considering cells of the same type 
and resistivity. The cell manufacturer was required by contract to use 
reasonable care in furnishing a simulator spectrum matching the Sun’s output 
so we assume minor spectral fluctuations had no effect. 

The a values used in the above equations axe calculated from the empirical 

formulae developed by Patterson et al (Reference 6). For 2 ohm-cm cells, 

2 x 2 x 0. 046 cm, I is expressed as 

sc 

I = C(T) H 
sc 
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where 


C(T) = c + c.T + mT 2 + c,T 3 + C/1 T 4 + c c T 5 
' o 1 2 3 4 5 


and 


c = 0. 914727 

o 


c 3 = 0. 22603-07 


c x = 0. 108713-02 


c 4 = 0. 171090-09 


c 2 = -0.695706-05 


c c = -0. 144039-11 

O 


(i3: 


<y is calculated from (13) by differentiating with respect to temperature. 
Representative families of such curves are shown in Figure 7. The value for 
a at 140 mW/cm and 28°C for these cells is 0. 11. 

Uncertainties in the voltage output of the solar cell are due to uncertain- 
ties in intensity, cell temperature and instrumentation accuracy. Voltage 
uncertainties can be calculated using the empirical Patterson formulae. The 
open circuit voltage dependence of a solar cell is expressed as 

V qc = A(T) + B(T ) logH 


where 


A(T ) 


a + a.T + a.T 2 + a.T 3 + a T 4 + a c T 5 
o 1 2 3 4 5 


16 


JPL Technical Memorandum 33-673 



and 



B(T) = b Q + bjT 

+ b 2 T 2 + b 3 T 3 

+ b 4 T 4 + b 5 T 

a o 

= 0.480654+03 

b o 

0.801129+02 

a l 

= -0.241702+01 

b l = 

0. 118039+00 

a 2 

= 0. 294965-02 

b 2 

-0. 224011-02 

a 3 

= -0.206159-04 

b 3 

0. 140478-04 

a 4 

= -0. 184560-06 

it 

& 

0. 816163-07 

a 5 

= 0.114636-08 

b 5 = 

-0. 548025-09 


The temperature coefficient of open circuit voltage is derived from (14) 
by differentiating: p = (9 V qc / 9T)^ an<3 tb e intensity dependence of voltage 

is given by 



B(T) 

H in 10 


These coefficients apply to 2 ohm-cm cells, 2x2 cm x 0. 46mm thick and 
strictly apply only to unirradiated cells. 

The V uncertainty is then 
oc J 



(15) 


where V qc is the open circuit voltage of the cell and dM is the uncertainty of the 
voltmeter calibration. 

Families of curves showing SV^^/OT and 9 V qc / 9H based on differentiating 
the Patterson formulae (14) are shown in Figures 8 and 9. 
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We use (15) to calculate an uncertainty dV Qc , then assume this voltage 
uncertainty applies over the entire I-V ’curve, but the percentage uncertainty 
varies over the curve as the inverse of the voltage. For dT = ± 1°C, 
dH = 4. 26 mW/cm^ (3. 04 percent of 140 mW/cm^) and dM = ± 1 mV (an 
uncertainty of ± 1 in the least significant digit of the digital voltmeter), 


dV 


oc 


= IT, 


245 x IT + (0. 254 x 4. 26) + (1) 


1/2 


mV 


= 2. 69 mV 


2 

for a measurement at 28°C and intensity of 140 mW/cm . Applying this voltage 
uncertainty to the 485 mV point, we compute a percentage uncertainty there of 


dV _ 
V 


0. 55% 


The uncertainty in voltage quoted here is a function of several factors 
which nearly all arise because of an inability to repeat measurements and 
experimental conditions exactly from cell to cell. A 0. 55 percent band exists 
on either side of the 485 mV setting. We should be willing to wager odds of 
19 to 1 that any particular cell had its voltage output within that band when its 
current output reading was recorded. Since the band is finite, and since the 
slope of the I-V curve is approximately -0. 25 mA/mV for 2 ohm-cm cells 
under these measurement conditions at the 485. mV point, this uncertainty 
in voltage will cause an additional uncertainty in the current reading of 
-0. 25 x 2.69 = 0. 67 mA or approximately 0. 56 percent. The nature of 
this uncertainty is to define an uncertainty area superimposed on any given 
point on the solar cell I-V curve. Such an area is illustrated in Figure 10 
for a portion of a typical 2 ohm-cm Mariner '71 cell near the 485 mV point. 
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To properly calculate a 95 percent confidence limit we should add the current 
uncertainties due to the direct current uncertainty expressed in equation (12) 
and the indirect current uncertainty arising from the voltage uncertainty 



- (3. 874) 2 + (0. 673) 2 
Al,p = 3. 93 mA 
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SECTION V 


UNCERTAINTIES IN THE SOLAR PANEL CALIBRATION 
AT TABLE MOUNTAIN 

When a panel has been assembled and has passed all its flight acceptance 
tests it is taken to Table Mountain, California, for the purpose of sun- 
light testing. It is instrumented to read panel temperature, pointed at 
the sun, and the entire I-Y curve for the panel is recorded. A balloon 
flight standard cell (BFS) is used during the test to establish the incident 
solar intensity during the test. The panel I-Y curve is used for categorizing 
the panels as to which will be flight panels (those best matched) and which will 
be flight spares. The panel I-V curve, which is the sum total output of an 
array of n cells in parallel by m cells in series, is compared with the I-V 
curve of a single solar cell of average output characteristics to calculate the 
manufacturing losses. Finally, the Table Mountain I-V curve is used to 
compute a standard I-V curve (presently at 145 mW/cm and 6 0°C) of a single 
cell on the panel by dividing the panel currents by n and the panel voltages by 
m. This standard curve is then used in the solar panel prediction program from 
which the panel behavior is computed for all expected solar intensities and 
panel temperatures during its space mission. The accuracy of the prediction 
program will be discussed in a later section. 

The uncertainties of the Table Mountain calibration arise from the 
original BFS calibration (± 1 percent), the temperature uncertainty of the BFS 
in use at Table Mountain (± 2°C), the panel temperature uncertainty (± 2. 5°C) 
and the distribution of temperature across the panel (unknown). Uncertainties 
in the calibration of the current and voltage measuring equipment will also 
contribute. The uncertainty in the intensity measurement is calculated by 
using equation (7), equation (2) and inserting an additional term for instrumenta- 
tion uncertainty. Assume the instrumentation uncertainty to be ± 1 digit 
in the least significant place or ± 1 mV out of 6 0 mV (0. 16 7 percent), and 
recall the ± 1°C in the BFS temperature contributed 0. 1 percent to compute: 


dH 

H 


[jo. 01) 2 + (0. 002 ) 2 + (0. 00167) 2 J X 100% 
1. 033% 
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In this case we can read the BFS cell with more assurance with the digital 
voltmeter because we do not have to contend with a dancing xenon arc as we do 
with the solar simulator. 

The uncertainty in the panel current may be computed by modifying 
equation (6) to account for the n cells in parallel on the panel or panel section 
and adding a term for instrumentation accuracy, dM/M. 



The voltage uncertainty is computed by modifying equation (15) to account for the 
m cells in series on the panel or panel section and adding the term for instru- 
mentation accuracy: 




m 

r-\ 

V 

oc 

3H ) 




(17) 


If we assume an instrumentation uncertainty of 0. 5 percent for both the 

voltage and current meters, and use appropriate test data taken from the 

July 17, 1970 Table Mountain testing of Mariner 1 7 1 panel S/N 019, Section E, 

(n=10, m=78, I =1.189 amperes, V = 41.66 volts, H = 127.3 mW/cm 2 , 

SC oc 

- 47. 2 C), we compute 


dl 

I 

P 


11 / 2 

(0. 01033) 2 + /L Q . x ... 0 ^ Q85 iL Z ,-- 5 \ + (0. 005) 2 
\ 1. 189 x 10 3 / 


X 100 


dl 

= 1. 161% (18) 
P 
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and 


dV 


oc 


V 


oc 


78 (-2. 28) 2 5\ 2 + / 78 (0. 277) 01033 x 12?> 3) ' 

, 41.66 X 1CT f \41.66 X 10 i 


— 1 1/2 

+ (0. 005) 2 X 1 00 = 1. 18% (19) 

2 

where we have used appropriate values for a, p, and dV id H at 127. 3 mW/cm 
and 47°C of 0.085 mA/°C, -2.28 mV/°C, and 0.277 mV/mW/cm respectively, 
which again, were obtained from figures 7, 8, and 9 or from differentiating the 
Patterson expressions. 

At the completion of the Table Mountain testing of a solar panel, the 
manufacturing losses are calculated. This is done by applying temperature 
and intensity corrections to correct the panel output to the temperature and 
intensity conditions at which the individual cells were measured. The panel 
I-V curve is transformed to a single cell curve by dividing the current by the 
number of cells in parallel and the voltage by the number of cells in series. 

The maximum power point of this transformed curve is compared with the 
maximum power point of the average single solar cell used in the panel. A 
percentage degradation is calculated and referred to as. the manufacturing loss. 
Spacecraft power people quote from 3 to 10 percent manufacturing loss. JPL 
typically finds a 6 percent loss. Several comments are in order. The panel 
current output is measured more accurately at Table Mountain than is the 
output of the single cells because the light intensity is better known. This may 
be mitigated by a less accurate solar spectrum at Table Mountain due to 
absorption in the atmosphere of selective portions of the spectrum. The voltage 
output is less uncertain for the single cells because the temperature is known 
more precisely. A part of the loss may be simply due to the uncertainties in 
intensity, temperature, and instrumentation. It is assumed here that the same 
BPS is not used at Table Mountain as is used at the Cell Vendor's establishment 
to calibrate his SWS cells. Neither are the same voltmeters used for measuring 
voltage and current at the two sites. Some of the uncertainties would be 
reduced if the same voltmeters and BFS cells could be used at each site. This 
procedure would not eliminate the absolute uncertainties in either measurement. 
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but would eliminate constant offset type errors in making the comparison, and 
would make the comparison between the single cell measurements and panel 
measurements more meaningful. 
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SECTION VI 


PREDICTING PANEL OUTPUT IN SPACE 

A. CORRECTION FOR SOLAR INTENSITY AND PANEL TEMPERATURE 

Predicting the panel output in space involves two types of corrections. 

The first is the correction for changing solar intensity and panel temperature. 
The second is for the effects of various environmental elements which damage 
the panel, namely electron and proton radiation damage and possibly UV 
degradation of the coverslide system. (Reference 7) 

The prediction of solar panel output in space is made using the JPL 

computer program Ml 32000. The calculational technique is based on the 

formulas of Sandstrom (References 8 and 9)> on some unpublished work by 

Blake (Reference 10), and miscellaneous adjustments by several others. The 

various temperature and intensity coefficients in the program are derived 

from the parametric data of Yasui (Reference 11). A reference I-V curve is 

stored in the program which is applicable to the single solar c ell/cover slide 

configuration of interest. Presently the I-V curve is for a Heliotek N/P 

2x2 cm cell, 0.46 mm thick, 2 ohm-cm resistivity, with a 7940 fused silica 

coverslide 0. 51 mm thick. The coverslide was mounted with RTV 602 adhesive. 

This is the configuration flown on Mariner *69 and *71. The reference I-V 

2 o 

curve is for a solar intensity of 145 mW/cm and a cell temperature of 60 C. 

It was derived from Mariner '71 panel measurements at Table Mountain 

(Reference 12) and therefore has all the inherent uncertainties of the Table 

Mountain measurements previously discussed as well as the uncertainties in 

correcting the experimental data to 145 mW/cm and 60 C from the Table 

Mountain data (taken at approximately 125 mW/cm and 50°C). Therefore all 

data on solar panel output from the program will have the uncertainties present 

that are discussed in this section, even any so-called uncorrected panel data 

2 

at 145 mW/cm and 6 0 C. The reference I-V curve consists of 30 I-V pairs 
spanning from well into the second current-voltage quadrant through the first 
quadrant and into the fourth quadrant. 
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The I-V curve for an arbitrary temperature intensity combination is 
found by translating the reference curve by the following equations: 

“(J) = I( J > + ^ref)(^ -l) + n* 2 (T 2 - Tl ) 

AI = II (J) -I(J) 

VI(J)= V ( J) -p 2 m (T 2 -Tj) - AI(R s2 + 0. 114) f ' ^ S (T 2 

( 20 ) 


Here 


1(J) and V(J) are current and voltage points on the reference I-V curve 


II (J) and VI (J) are current and voltage points on the transformed 
I-V curve. 

I (ref) = reference curve short circuit current 

S C 




H, 


R 


s2 


new solar intensity 
new cell temperature 
reference curve solar intensity 
reference curve cell temperature 

averaged short circuit temperature coefficient at and 
referenced to 

averaged open circuit voltage temperature coefficient at 
and referenced to 

solar cell series resistance at T^ 

solar cell curve shape correction factor at intensity H_ 

c* 
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m = number of cells in series on the panel 
n = number of cells in parallel on the panel 

These equations simply translate the whole I-V curve parallel to the 
current axis by the same amount that the short circuit current is translated 
with varying temperature and intensity. The I-V curve is translated parallel 
to the voltage axis by the same amount that the open circuit voltage is trans- 
lated with temperature and intensity and by an additional amount equal to the 
1R drop across the panel calculated for the new current. If the equation for a 
single cell is desired, m and n are set equal to 1 and the factor of 0. 114 is not 
added to R 7 , The 0, 114 appears to be an adjustment found necessary to 
account for interconnect and panel harness resistance, and should probably 
be temperature -dependent. 

Note that the "temperature coefficients" a and (3 are derived relative to 
60°C only (reference 10), and must be changed if a reference curve at a 
different temperature is desired. These coefficients in their present form 
have been calculated in an approximate fashion and should be updated. The 
cell series resistance derivation was done by a standard technique originally 
suggested in reference 13. The technique must be carried out very carefully 
to avoid certain pitfalls and snares. The resistance is calculated using points 
from two I-V curves of the same solar cell illuminated at different intensities. 
The calculated resistance varies according to what illumination intensities 
happened to have been chosen by the experimenter. Typically the cell resist- 
ance measured in this manner will vary by about ±0. 05 ohm for the same cell, 
depending on the intensity levels used. The measurement should properly be 
done by examining a large number of cells and a least squares fit made to all 
the data in order to derive a temperature dependence and a 95 percent confi- 
dence envelope. The accuracy of the cell resistance values in use is not known 
nor are the confidence limits known. 

After making temperature /intensity corrections using the a's, p's, and 
R ? , it was found that the knee of the corrected I-V curve was improperly 
rounded. The curve shape correction coefficient, K^, was therefore included 
in formula (13) and was empirically adjusted to give the "best fit" (Reference 10) 
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to the laboratory parametric I-V curves. The criterion for establishing best 
fit is not stated. Neither are confidence limits established. The factors 
were found to be mainly a function of solar intensity. 

The derived a's, p's, R f s and K 's were used to recalculate (Refer- 

S Li Ca 

ence 10) the original parametric data. Here transformations were made from 
a reference curve at 253 mW/cm solar intensity and 60°C. A table of cor- 
relations was calculated. The correlation table appears approximately correct. 

Maximum deviations of the calculated I , V , and P from the laboratory 

sc oc max 7 

values were 1.36 percent, 1.36 percent, and 1.63 percent respectively, but 
these values may be optimistic. The optimism may be expressed by using the 
maximum deviations as 95 percent confidence limits. The original para- 
metric data of Yasui were obtained with a sample size of 13 solar cells. He 
computed statistical parameters for the cells including 95 percent confidence 
limits. Considering only the temperature - solar intensity regime bounded by 
-20°C, 60°C, 50 mW /cm^ and 250 mW/cm^ which is of primary interest 

here, maximum 95 percent confidence limits of 0,97 percent and 0.33 per- 
cent were found for I and V respectively. Prediction errors found by 
Blake over this same temperature -solar intensity regime are 1. 13 percent for 
I sc an d 0.78 percent for V . The uncertainty of the predicted solar panel 
current may be expressed as 


dl 

I 

P 


dH 

H 


a dT 
o t 


nndT 


TM 


4* S T + e. 


( 21 ) 


where S. j. is the statistical 95 percent confidence limit of the parametric current 
data and is the prediction error in current due to the approximate nature of 
the computing technique. dT is the temperature uncertainty of the cells in the 
test chamber during the parametric study (± 1°C) with a corresponding a and 
current 1^., dT_^ is the solar panel temperature uncertainty with an a appropriate 
to (H, T }. n is the number of cells connected in parallel on the panel or panel 

ir 

section. (dl/I)^^. is the uncertainty in measuring the solar panel current at 
Table Mountain. This term must be added in if the reference I-V curve in 
the computer program is derived from the Table Mountain data as is the 
present practice. 
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The solar intensity uncertainty in equation (21) consists of two parts. 

One is the uncertainty of the simulator intensity when generating the parametric 
data and is calculated by using appropriate values in equations (8) and (9), 
namely 1.0 percent for the BFS calibration uncertainty, 1.0 percent for the 
uncertainty due to simulator drift, 0. 1 percent for the uncertainty due to BFS 
temperature control and measurement, and 0.5 percent for the uncertainty 
associated with the calibration and reading of the digital voltmeter. The total 
intensity uncertainty of the simulator when generating parametric data there- 
fore is: 



(0. 01) 2 + (0. 01) 2 + (0. 001) 2 + (0. 005 ) 2 


1/2 

x 100 = 1. 50% 


The other part comes from the uncertainty in solar intensity due to uncertain- 
ties of spacecraft position with respect to the sun and the incidence angle of 
sunlight on the panels. Uncertainties in the solar constant do not come into 
play here because all measurements are ratioed to balloon flight standard 
cells which were measured at air mass zero at one solar constant, whatever 
its value. The solar intensity incident on the panel is calculated by 


K cos 8 

2 

r 


( 22 ) 


where r is the spacecraft to sun distance, 6 is the incidence angle, and K is a 
constant. 



2 H t 

- tan 0 d0 

r 


(23) 
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On Mariner '71 d0 was about 1/4 degree (4. 36 x 10"^ radian) and dr was 

g 

approximately 10 km. r is of the order of 10 km. Therefore 



r ' \Z 
2 X 10) 

io s ) 


-|l/2 


-3 


-3.2 


+ (4. 36 X 10 X 4. 36 X 10 ") 


X 100 = 0, 0019% 


and the combined intensity uncertainties in solar intensity are 


dH 

H 



1. 50% 


(24) 


It is apparent that the uncertainty in the solar intensity incident on the space- 
craft solar panels is negligible in comparison to the other uncertainties. 

As an example of the application of equation (2l), assume we want to 

2 

find the uncertainty in the predicted solar panel current at 100 mW/cm solar 

intensity and panel temperature of 0°C, The appropriate a and I for these 

s c 

conditions are 0. 11 mA/ C and 90 mA (for a single 2 ohm- cm solar cell). 
Assume a solar panel temperature uncertainty of 5°C and use the value 
computed in equation (18) for the relative current uncertainty in the Table 
Mountain measurements (0. 01161) and calculate a total uncertainty in the 
predicted current as follows: 


dl 

_2 = 
I 

P 


(0. 015) 2 + '°* 11 X 1 


90 


r * ( ? 


. 1 1 x 5 

90 


+ (0. 01 16 1 ) 2 + (0. 0097) 2 


|l/2 


+ (0. 0113)' 


x 100 = 2. 49% 


(25) 
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In a similar fashion, the uncertainty in voltage is derived from a 
modification of equation (15) to give 


dV 

oc 

V 


oc 




3V 

oc 

3H 



f av 


ocp 


ocp, 


+ s v + 


TM 


( 26 ) 


Sy is the 95 percent confidence limit for the parametric voltage data, e^. is 
the computer voltage prediction error, V q is the panel open circuit voltage, 
and 



is the uncertainty in measuring the solar panel voltage at Table Mountain. 

This term must be added in if the reference I-V curve in the computer program 
is derived from the Table Mountain data as it is in present practice. dH is 
the value calculated using equation (24), 

As an example of the application of equation (26), we will calculate the 

2 

uncertainty in the predicted solar panel voltage <at 100 mW/cm solar intensity 
and panel temperature of 0°C. Appropriate values for the required parameters 
are: V qc = 641 mV, dH = 0. 015 x 100 mW/cm 2 = 1.5 mW/cm 2 , (3V /3T) = 

-2. 18 mV / °C, dT t = 1°C, dT p = 5°C, (3 V qc / 3H) = 0. 35 mV-cm 2 /mW 

Sy = 0. 0033, ey = 0. 78 percent (Blake's "maximum" prediction error of V ), 

and (dV /V ) = 0. 0118 from equation (19). 

OCp OCp 


dV 


oc 


oc 


1 - 2 . 18 x l\ 2 , 1 - 2 . 18 x 5\ 2 ( 0 . 35 x 1. 5^ 2 , 

r 641 ) + v 6 4 1 ' / + \ 64l ) + (0.0118) 


+ (0. 0033) 2 + (0. 0078) 2 


1/2 


x 100 = 2. 26% 


(27) 
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Since the entire I-V curve is translated along both axes by the computer 

program and the amount of the translation has been computed using I and 

V behavior as the basis, errors in 1 and V apply to each point on the 
oc sc oc 1 r 

curve. An error of 1 mA in I , for example, is an error of 1 mA all along 

s c 

the I-V curve, and it is a higher percentage at all other points. Similarly 
for V . An uncertainty in power is computed simply as 



/dV 

V \ 

/dl 

1 V 

/ oc 

oc ] + 

< sc 

sc \ 

\ V 

v ) I 

l I 

I / 

\ oc 

‘■nP/ 

\ sc 

mp/ 


(28) 


at the maximum power point. Note that this is in effect an approximation 
of the four sided area of Figure 10 by a rectangle. 

The hypothetical Mariner *71 panel section for which we have evaluated 

formulas (21) and (26) at 100 mW/cm and 0 C (an m x n array of 78 x 10 cells) 

would have a maximum power of 35. 5 watts at 42. 26 volts and 0. 839 amperes 

under these conditions computed by using the parametric data of reference 11 

but ignoring manufacturing losses. Correspondingly, approximate I and 

s c 

V values would be 0. 92 amps and 5 0, 0 volts, dl = 0. 92 x 0, 0249 = 
oc sc 

0. 0229 amps and dV = 50. 0 x 0. 0226 = 1. 13 volts. Applying equation (28), the 
uncertainty in maximum power due to the prediction process is 


dP 


m 


m 



x 100 = 3. 82% 


(29) 


This is a typical power prediction uncertainty for Mariner type solar 
panels operating between Venus and Mars when the panel components have been 
measured using today's techniques. A summary of the contributions to the 
error expressed in equation (28) is shown in Table I. 

Squaring and adding the values in the columns to the right will give 
relative uncertainties in 1^ and V respectively. The examples are for a 
calculation of panel performance at 100 mW/cm and 0°C. Other numbers 
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Table I. Sources of Solar Panel Uncertainty 


Parameter 

Uncertainty 






I. I 

Sc 

SHORT CIRCUIT CURRENT 




dH 

(Simulator) 

BFS calibration 
BFS temperature 
DVM calibration and read 
Simulator drift 

0. 015 
0. 01 
0. 001 
0. 005 
0. 01 



dH 

(Solar Intensity) 
S/C - Sun distance 
Pointing error 

1.9 - 05 

2.0 - 07 

1.9 - 05 



dH 

(Total) 



0. 015 

dT t 

(Parametric Data)* 

1. 0°C 




/edT A 

1 y 1 } at < 100 > °) 



1. 22-03 

dT 

P 

(S/C panel temperature)* 

5. 0°C 




/odT \ 

\ -p- 2 ] at (100, 0) 



6. 11-03 

dl/I 

(Table Mountain Reference Curve) 



0. 01161 

S I 

(Parametric statistics) 



0. 0097 

e I 

(M132000 approximation) 



0. 0136 

II. V 

oc 

OPEN CIRCUIT VOLTAGE 




dH 

(Simulator + Solar Intensity) 





/ aV oc\ dH 
\ 3H / V qc at (100, 0) 

0. 015 


8. 19-04 

dT t 

(Parametric Data)* 

1°C 




/SV \ dT 

— t 

\ 9T /„ V at (100, 0) 



3.40-03 

dT 

P 

(S/C panel temperature)* 

5°C 




/av \ dT 

iTf) V ~ ’ 1 at (100, 0) 
H oc 



0. 017 

dV 

oc 

V 

oc 

(Table Mountain Reference) 



0. 0118 

s v 

(Parametric statistics) 



3. 30-03 

e v 

(M1320QQ approximation) 



1. 36-02 


Starred data ia absolute, other non-starred data is relative (xxx, yy) means at solar 
intensity of xxx mW /cm , and yy°C. The quantities in the right column are based o: 
derived from the quantities in the middle column. 
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must be substituted for the five values in the table written with (100, 0). The 
values for e, and e in the table are pertinent to the temperature- solar 
intensity regime -40 to 140 C, and 25 mW/cm to 540 mW/cm . Lesser 
values than these were used in the example calculated because we were there 
restricting ourselves to a smaller temperature- solar intensity regime, and the 
maximum value for the restricted regime was used. After finding the relative 
uncertainties in panel voltage and current, equation (28) is used to find the 
uncertainty in maximum power. 

B. CALCULATION OF A MISSION POWER PROFILE 

This section will treat the calculation of solar panel output during the 
course of a mission, allowing for various environmental degradation factors 
and the changing temperature and solar intensity. An initial I-V curve for the 
panel just after launch is assumed. If the panel is to be manufactured from the 
same kind of solar cells as was a previously built well known panel, the start- 
ing I-V curve may be taken from that panel. This initial curve may be 
obtained from the I-V curve of a single solar cell with the voltage multiplied 
by m, the current multiplied by n and the resultant curve modified by some 
sort of panel assembly loss factor. After a panel has been built and tested at 
Table Mountain a better estimate of an initial I-V curve may be obtained by 
correcting to initial space conditions using the techniques of Ml 32000. This 
initial I-V curve thereafter is translated parallel to the voltage and current axes 
by the program to take into account the known cell response to the changing envi- 
ronments. This procedure has the inherent assumption that the environments 
produce little change in the temperature and intensity coefficients of voltage and 
current. It also assumes there is no significant fluence of low energy protons 

incident on unprotected cell surface areas. (Significant fluences would be, for 
112 

example, 10 p/cm of 0.2 MeV protons incident a 10 mil gap near the bus bar. ) 

The first step in calculating the panel output is to calculate the response 

curve as a function of mission time assuming no environmental degradation. 

This assumes that the time variation of the panel temperature and the incident 

solar intensity is known. The uncertainties in I , V , and P may be 

sc oc max 3 

assigned to the calculated time dependent curves at whatever confidence level 
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desired by using the principles discussed in this paper. The result of the 
Ml 32 000 calculation may be considered to be the predicted average I-V curve 
as a function of time and the uncertainty calculation will establish a band of 
curves above and below the average curve, each band corresponding to a 
particular confidence limit. 

The non-degraded profile curve will in general require modification to 
account for degrading environments. A degrading environment is considered to 
be one which produces an irreversible and usually deleterious effect. The 
basic technique is to consider each environment individually and calculate its 
modification to the current and voltage of the panel. For example, suppose 
the panel must fly through two environments, A and B, during the mission. 

Also suppose environment A alone would reduce the solar cell short circuit cur- 
rent by percent, and the open circuit voltage by percent, and environ- 
ment B alone would reduce the short circuit current by Xg percent and the open 
circuit voltage by percent. Then the resulting I-Y curve can oe found by 
translating the undegraded I-V curve. Translate the current by an amount 



and the voltage by an amount 



(31) 


to get The new I-V curve. I and V are the undegraded short circuit current 
and open circuit voltage of the panel assumed to be corrected for intensity and 
temperature. The maximum power will be determined from the resultant curve. 
If more than two degrading environments are involved, additional multiplicative 
terms are appended to equations (30) and (31). 

The degrading environments likely to require attention are radiation and 
the mysterious monotonic degradation (her einafter referred to as MM) which is 
seen on some spacecraft and is apparently independent of bulk damaging radia- 
tion (Refs. 7 and 16). Other environments that may need attention on some 
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missions are micrometeoroids, thermal cycling, and ultraviolet interaction 
with coverslide adhesives or optical coatings. 

For satellites orbiting the Earth in the Van Allen belts, the radiation flu- 

ence increases more or less monotonically as the spacecraft stays in Earth 

orbit with an occasional fluence of solar flare protons added from time to time. 

The interplanetary missions flown by JPL so far have not been subjected to 

radiation fluences of consequence. These spacecraft receive only minimal 

radiation doses as they pass through the Van Allen belts. Their major 

radiation exposure results from the sporadic emanations from the Sun when it 

emits protons during a solar flare. Therefore the JPL panel designer does 

not usually have a smooth curve detailing radiation exposure as a function of 

time. Instead, he will most likely have a curve which plots proton fluence as 

a function of the probability of receiving less than that fluence. A fluence- 

energy spectrum definition will also be given in some form. It is usually up to 

the individual panel designer to select the radiation risk he wishes to include 

in his design, by selecting a probability level on the fluence curves. He must 

also assume that the total radiation dose will occur all at once, since nearly 

the entire exposure may be due to the protons emitted during a single solar 

event. He must further assume that this event will occur at the least optimum 

time, probably when the spacecraft is closest to the sun since the fluence goes 

2 , 

roughly as 1/r , or perhaps just prior to or during the peak power demand 
of the spacecraft. 

To properly factor the radiation damage into the panel power prediction, 
the radiation environment for 5 0 percent probability of occurrence should be 
selected. This amount of radiation should be coupled with the standard 
techniques for computing solar cell degradation (References 14 and 15) and a 
percentage remaining P computed using equations (3 0) and (31). This 
fraction should then be multiplied by the predicted time dependent performance 
curve for a non -irradiated panel. The 95 percent confidence limit of panel 
degradation due to radiation may be established as follows: Take the radiation 
level on the fluence vs probability curve at the point where the probability is 
95 percent that the spacecraft will receive no more than this level. Calculate 
the cell degradation and the percentage electrical power remaining as before. 
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Multiply this fraction by the predicted time dependent performance curve for a 

non- irradiated panel. The difference in power defined by subtracting the 

powers calculated using the twofluences gives a dP f which may be expressed 

as a relative power, squared, added to the square of the power uncertainty 

calculated in equation (28), and square rooted to give an overall uncertainty. 

Similar procedures may be used for I and V and for other confidence 
r } sc oc 

limits. 

This procedure should be followed for factoring in the panel degradation 
for other environmental agents. Establish an average degradation due to the 
effect of the environment, then find some way of calculating a degradation band 
encompassing the desired confidence level. If the band width is expressed in 
maximum power, then the uncertainty may be added to the previously computed 
maximum power uncertainty by computing the relative uncertainty (using half 
the band width), squaring, and adding to the square of the uncertainty due to 
other effects and taking the square root of the total to give the total uncertainty. 

If the designer has used an accurate starting I- V curve and followed the 
above steps he has now effectively arrived at a family of curves which give 
solar panel output as a function of mission time. Each curve of the family 
corresponds to a different level of confidence that the panel will, meet or exceed 
the plotted power level. The whole family of curves may be moved up or down 
(translated parallel to the power axis), by simply adding or subtracting the 
total number of cells on the panel to make slight adjustments in the power 
profile. Large translations will not be accurate because of the differing ways 
in which m and n affect the power uncertainty of equation (28). The family of 
curves may be superimposed on the required power profile curve for the 
mission and translated up or down until the panel power output curve exceeds the 
power requirement curve at all times in the mission. Comparison w r ill be 
made at the confidence level dictated by the project philosophy. If a substantial 
translation of the original calculated curve has been necessary to give the 
proper sizing, a recalculation of the family of curves will be necessary with 
the new values of n and m. Of course, if there is only one point in the mission 
where the power demand and power available curves may cross, then the 
entire set of calculations required to size the panel may be performed at that 
point. Having established n and m, the calculation of the family of power 
profile curves may be completed. 
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The mysterious monotonic radiation has been observed on the IDCSC 
spacecraft flying at synchronous altitude (Reference 7) and on Mariners 6, 7, 
and 9 (Reference 16). Patterson reports in reference 16 that the MM degrada- 
tion on these Mariners has been fairly repeatable. The MM degradation has 
been referred to as uY degradation, mainly because better explanations have 
not been forthcoming. Other possibilities exist. All the IDCSC spacecraft 
used a primer and the Dow XR-6-3489 adhesive in mounting the coverslides. 
RTV-602 adhesive was used to mount the Mariner coverslides, but no primer 
was used. The IDCSC spacecraft used 10 ohm- cm Heliotek cells, the Mariners 
used 2 ohm- cm Heliotek cells. There is enough difference in the two solar 
cell systems so that the MM degradation seen on the two spacecraft could arise 
from different causes, or a combination of the same causes but a different mix. 
For example the IDCSC cells may be degrading because of a true uV degradation 
in the cover slide primer. Both spacecraft may suffer additional degradation 
due to uY degradation of the adhesive or optical coatings. Both may very 
likely be degrading because of contamination on the surfaces of the coverslides 
which darkens in the presence of radiation. Contaminants may arise from the 
spacecraft itself or from the vacuum chambers where the spacecraft were 
environmentally tested. An Air Force Tactical Communications Satellite 
(TACSAT) built by Hughes Aircraft Company did not experience the MM degrada- 
tion. On this spacecraft, all the spaces between solar cells are fully grouted 
with an epoxy, and this panel was thoroughly cleaned with a mild abrasive 
after environmental testing and prior to launch. Surface contamination from 
the environmental vacuum chamber was probably removed quite effectively, 
but the cleaning was severe enough so that the uV cutoff filters on the cover- 
slide could have been injured or removed. (Reference 17) 

At any rate the MM degradation which has been repeatable, may not be 
repeatable or predictable in the future if panel assembly and environmental 
testing procedures are changed. The cause of this degradation should be 
investigated. 
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SECTION VII 


RECOMMENDATIONS 

1. Coverglasses should be applied to solar cells prior to the series of solar 
cell measurements to be used for matching the cells into submodules. 

2. An investigation should be launched to determine the reason why there is 
so much scatter in the solar cell short circuit current loss due to glass- 
ing, and to try to find methods of reducing the scatter. One possible 

method is the use of TiO antireflection coatings on the solar cells instead 

x 

of SiO. 

3. Rewrite the JPL Solar Panel Prediction program (Ml 32000). The 

revised program should be tied more directly to the original parametric 

data generated in the simulation laboratory. For example it would he 

better to have a table of parametric data for I as a function of intensity 

s c 

and temperature from which I at all intermediate intensities and 

sc 

temperatures could be interpolated. As it is now, several factors must 
be calculated from the parametric data and then put into Ml 32000 before 
the program can be used. The program should have the capability of 
calling a radiation degradation subroutine to calculate the effect of any 
radiation environment on the electrical parameters of the panel. Other 
subroutines may be added, for example, a routine to compute cover slide 
darkening as a function of uY exposure if it can be proven that a functional 
relationship exists. 

4. Tighten up the specification on setting up the simulator intensity at the 
solar cell vendor's establishment. The procedures currently in use 
give an enormous uncertainty in the simulator intensity. 

5. Investigate carefully the exact cause of the mysterious monotonic degra- 
dation of the solar panels and the I - I transducers on the Mariner 

s c s c r 

solar panels. 
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Fig, 1. Current at 0,485 volts, distribution for 2-ohm-cm cells 



124 126 128 130 132 134 136 138 140 142 144 


*$c r mA 

Fig. 2. Short circuit current distribution of 10-ohm-cm solar cells 
before and after glassing with 6 -mil micro sheet 
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LOSS DUE TO APPLICATION OF COVERSL1DE, mA 
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^sc, mA 

Fig* 3. Short circuit current distribution of 2-ohm-cm solar cells 
before and after glassing with 6 -mil 7940 fused silica 
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INITIAL 5H0RT CIRCUIT CURRENT, mA 

Fig. 5. Current reduction after cover slide mounting (2x4 cm cells) 



Fig. 6. 95% current confidence band on a 2-ohm -cm solar cell 
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Fig. 7. a = (3 I gC / 9 T)^; based on differentiating Patterson's expression 
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Fig. 8. (3V oc /9T)pj- based on differentiating Patterson's expression 
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Fig. 9. 
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based on differentiating Patterson’s expression 
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Fig. 10. 


Portion of 2-oh.m-cm solar cell I~V curve near P 
for AMO, 28*C 
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